Table 2: The numerical values of coefficients entering the expression of as defined in (4.13). 
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Abstract. This paper investigates the non-linear self-interaction of quadrupole 
gravitational waves generated by an isolated system. The vacuum Einstein field 
equations are integrated in the region exterior to the system by means of a post- 
Minkowskian algorithm. Specializing in the quadrupole-quadrupole interaction (at 
the quadratic non-linear order), we recover the known results concerning the non- 
local modification of the ADM mass-energy of the system accounting for the emission 
of quadrupole waves, and the non-local memory effect due to the re-radiation by 
the stress-energy distribution of linear waves. Then we compute all the local 
(instantaneous) terms which are associated in the quadrupole-quadrupole metric 
with the latter non-local effects. Expanding the metric at large distances from the 
system, we obtain the corresponding radiation-field observables, including all non-local 
and transient contributions. This permits notably the completion of the observable 
quadrupole moment at the 5/2 post- Newtonian order. 
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1. Introduction 

In general relativity, the multipole moments of any finite distribution of energy 
and momentum interact with each other in vacuum, through the non-linearities of the 
field equations. In particular, the multipole moments which describe the gravitational 
waves emitted by an isolated system do not evolve independently, but rather couple 
together (including with themselves), giving rise to non-linear physical effects. 

The simplest multipole interaction which contributes to the radiation field is that 
between the (mass- type) quadrupole moment and the mass monopole M. The 
latter moment is the constant mass-energy of the source as measured at spatial infinity 
(ADM mass). Associated with the multipole interaction M^j x M is the non-linear 
effect of tails. This effect is due to the backscatter of linear waves (described by 
My) onto the space-time curvature generated by the mass-energy M. The tails can 
be computed within the theory of gravitational perturbations of the Schwarzschild 
background (see e.g. [1-3]). A consequence of the existence of tails is the non-locality 
in time, as the tails are in the form of integrals depending on the history of the source 
from — oo in the past to the retarded time t — r/c [4-7]. It is known that the tails 
appear both in the radiation field and in the radiation reaction forces at the 3/2 post- 
Newtonian order (1.5PN, or order c -3 when c — > oo) relatively to the quadrupole 
radiation [8,9]. 

Next in complexity is the interaction of the quadrupole moment with itself, or 
quadrupole self-interaction x Mj-i- Two closely related non-local (or hereditary) 
effects are known with this particular interaction. As shown by Bonnor and 
collaborators [10,11,4,6], the total mass M gets modified by a non-local integral 
accounting for the energy which is radiated by the quadrupole waves (there is 
agreement with the Einstein quadrupole formula). On the other hand the radiation 
field involves a non-local contribution whose physical origin is the re-emission of waves 
by the linear waves [12-16,9]. This contribution can be easily computed by using 
as the source of waves in the right side of the Einstein field equations the effective 
stress-energy tensor of gravitational waves (averaged over several wavelengths). As 
shown by Christodoulou [14] and Thorne [15] this implies a permanent change in the 
wave amplitude from before to afterward a burst of gravitational waves, which can 
be interpreted as the contribution of gravitons in the known formulas for the linear 
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memory [17,18]. The latter non-linear memory integral appears at the 2.5PN order in 
the radiation field (1PN order relatively to the tail integral). 

Now the metric which corresponds to the quadrupole-quadrupole interaction 
involves also, besides the non-local contributions, many terms which by contrast 
depend on the multipole moments at the sole retarded instant t — r/c. In the following 
we shall often qualify these local terms as instantaneous (following the terminology 
of [8]). The instantaneous terms in the radiation field (to order 1/r) are transient in 
the sense that they return to zero after the passage of a burst of gravitational waves. 
On physical grounds, it can be argued that the instantaneous terms do not play a 
very important role. However, these terms do exist and form an integral part of the 
field generated by very relativistic sources like inspiraling compact binaries (see e.g. 
[19]). In fact, the complete quadrupole-quadrupole metric including all the hereditary 
and instantaneous contributions will be needed in the construction of very accurate 
theoretical wave forms to be used by the future detectors LIGO and VIRGO. 

The present paper is devoted to the computation of the quadrupole-quadrupole 
metric, using the so-called multipolar-post-Minkowskian method proposed by 
Blanchet and Damour [20,21] after previous work by Bonnor [10] (his double-series 
approximation method) and Thorne [22] (how to start the post-Minkowskian iteration 
using STF multipole moments) . The instantaneous quadrupole-quadrupole terms have 
never been computed within the multipolar-post-Minkowskian method. However, they 
have been computed by Hunter and Rotenberg [6] within the double-series method, 
in the case where the source is axi-symmetric. The non-linear interaction between 
quadrupoles has also received attention more recently within the double-series method 
[23]. 

The main result of this paper (having in mind the application to astrophysical 
sources to be detected by VIRGO and LIGO) is the completion of the observable 
quadrupole moment of a general isolated source at the 2.5PN order. To reach this 
result we also take into account previous results concerning the tails at the 1.5PN order 
[9], and the multipole moments given by explicit integrals over the source to 2.5PN 
order [24,25]. (Note that the 2.5PN approximation in the observable quadrupole 
moment gives no contribution to the phase evolution of inspiraling compact binaries 
[25], however it will be required when we compute the 2.5PN wave form.) 

In a following paper [26], that we shall refer to as paper II, we shall investigate 
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the monopole-monopole-quadrupole interaction (at the cubic non-linear order) which 
enters the radiation field at the 3PN approximation. The present paper and paper 
II are part of the program of computing the field generated by inspiraling compact 
binaries with 3PN and even 3.5PN accuracy (see [27-30] for why such a very high 
accuracy is necessary). 

The plan of the paper is as follows. In Section 2 we summarize from [20,21] 
the method for computing the field non-linearities. In Section 3 we investigate 
(following [8]) the general structure of the quadratic non-linearities. Section 4 deals 
with the explicit computation of the quadrupole-quadrupole metric. The results are 
presented in the form of tables of numerical coefficients. Finally, in Section 5, we 
expand the metric at infinity and obtain the observable moments in the radiation 
field (essentially the quadrupole). The technical formulas for integrating the wave 
equation are relegated to Appendix A, and some needed results concerning the dipole- 
quadrupole interaction are derived in Appendix B. Henceforth we pose c = 1, except 
when we discuss the post- Newtonian order at the end of Section 5. 



2. Non-linearities in the external field 

We summarize the method set up in [20] for the computation of the quadratic and 
higher non-linearities in the field generated by an isolated system. The computation 
is performed in the exterior (vacuum) weak-field region of the system, where the 
components of the gravitational field W are numerically small as compared to one. 
Here h^ v denotes the metric deviation W = \f—gg* iv — t]^, with g^ u the inverse and 
g the determinant of the metric g^, and the Minkowski metric diag(— 1, 1, 1, 1). 
The field h^ v admits in the exterior weak-field region a post-Minkowskian expansion, 

h? v = Gh^ + G 2 h^ + ... + G n h^ + ... , (2.1) 

where G is the Newtonian gravitational constant, which plays here the role of book- 
keeping parameter in the non-linearity expansion. The first term in (2.1) satisfies the 
vacuum Einstein equations linearized around the Minkowski metric. In harmonic (or 
De Donder) coordinates this means two equations, 
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Uh^ = , 
d v h»" = . 



(2.2a) 
(2.26) 



In the first equation □ denotes the flat space-time wave (d'Alembertian) operator. 
The second equation ( diver genceless of the field) is the harmonic gauge condition. The 
equations (2.2), supplemented by the condition of retarded potentials, are solved by 
means of a multipolar expansion. It is known that only two sets of multipole moments, 
the mass-type moments Ml and current-type moments Sl (both depending on the 
retarded time t — r), are sufficient to parametrize the general multipole expansion (see 
e.g. [22]). In terms of these moments the general solution reads [22] 



£>0 



C = "4$:^/a i [r- 1 M L (t-r)] , 

r-'Mgl^t-r) 



(2.3a) 



£>\ 



i-yi 



+ 4 E77XTw £ -^-i[ r " 1 ^-i( t - r )] > 
e>i { + >' 



(2.36) 



£>2 



a 



L-2 



e>2 



r 1 e a b{ l S$ L _ 2 {t-r) 



(2.3c) 



Here the superscript (n) denotes n time derivatives. The index L is a shorthand 
for a multi-index composed of i indices, L = (similarly, L — 1 = i\i2---ii-\, 

aL — 2 = ai\...ig-2, and so on), and &l denotes a product of t space derivatives, 
9l = di x di^...di t . The multipole moments Ml and Sl are symmetric and tracefree 
(STF) with respect to all their indices (see [31] for a recapitulation of our notation 
and conventions). 

A priori the multipole moments M^(t) and SL(t) are arbitrary functions of time. 
They describe potentially the physics of a general isolated source as seen in its exterior 
field [22,20]. The only physical restriction is that the mass monopole M (total mass- 
energy or ADM mass), the mass dipole M; (position of the center of mass times the 
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mass), and the current dipole Si (total angular momentum) are constant. Technically 
speaking, this is a consequence of the harmonic gauge condition (2.2b) (see e.g. 
[22]). In this paper we shall set the mass dipole Mj to zero by shifting the origin 
of coordinates to the center of mass. In order to describe an isolated system, we must 
implement a condition of no-incoming radiation, ensuring that the radiation field is 
entirely generated by the system. We assume that the field is stationary in the remote 
past, i.e. that the moments Mi,(t) and Sx(£) are constant before some finite instant 
in the past, say when t < —T (hence Mj cannot be a linear function of time and is 
necessarily constant or zero). This assumption may seem to be somewhat restrictive, 
but we can check a posteriori that the formulas derived in this paper and paper 
II admit a well-defined limit when — T — > — oo in more general physical situations, 
such as the formation of the system by initial gravitational scattering. The multipole 
moments M L (t) and -Sx(t), subject to the previous restrictions, play the role of "seed" 
moments for the construction of the exterior field (2.1). In particular we shall express 
the results of this paper in terms of products of My with itself. We shall not use the 
expressions of the multipole moments Ml and Si, as explicit integrals over the source. 
However, these expressions are known in the post-Newtonian approximation [24,25], 
and should be used in applications. 

The coefficient of G 2 in (2.1) is the quadratically non-linear metric, whose precise 
definition we recall. The field equations for this coefficient read, still using harmonic 
coordinates, 



Uh% v = N% v , 
d„hii v = . 



(2.4a) 
(2.46) 



The d'Alembertian equation involves a quadratic source N!£ v = N^ v {h\, hi) generated 
by the linearized gravitational field (2.3), where 



N^(h, h)=- h pa d p d a h^ + -d tl h pa d u h p ' J - -d»hd u h 



- 2d^h pa d p h^ a + d a W p {d a h v p + d p h va ) 



- l -d x h pa d x h p ° + l -d p hd p h + l -d p h aX d°h pX 



(2.5) 
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From (2.4b) we deduce 



. 



(2.6) 



Then the quadratic metric h^" , solving (2.4) and the condition of stationarity in the 
past, is obtained as the sum of two distinct contributions, 



Basically the first contribution u^" is the retarded integral of the source N$ v '. 
However, in our case the source is in the form of a multipole expansion (valid only 
in the exterior of the system and singular at r = 0), so we cannot apply directly the 
usual retarded integral operator, whose range of integration intersects the system at 
retarded time. A way out of this problem, proposed in [20], consists of multiplying 
the actual source term by a factor (r/ro) s , where B is a complex number and 
r denotes a certain constant having the dimension of a length. When the real part 
of B is large enough all the power-like (except for logarithms) singularities of the 
multipole expansion at the spatial origin of the coordinates r = are cancelled. 
[Actually we consider separately each multipolar pieces, with given multipolarities 
£, so that the maximal power of the singularities is finite, and B can indeed be chosen 
in such a way] Applying the retarded integral on each multipolar pieces of the product 
{f /fo) B ' results in a function of B whose definition can be analytically continued 
to a neighbourhood of B — 0, at which value it admits a Laurent expansion. The finite 
part at B = (in short FPs =0 ) of the latter expansion is our looked-for solution, as 
it satisfies the correct wave equation (Du^ = N% v ), and is like in the form of 
a multipole expansion. Note that the latter process represents simply a convenient 
mean to find a solution of the wave equation whose source is in the form of a multipole 
expansion. Other processes could be used as well, but this one is particularly powerful 
as it yields many explicit formulas to be used in practical computations (see Appendix 
A below and Appendix A of paper II). Hence the first contribution in (2.7) reads 




(2.7) 




(2.8) 



where U R denotes the usual retarded integral 
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(□-V)(x,t) = -i- HI ^^/(x',t-|x-x'|). (2.9) 

When dealing with the metric at quadratic order, it can be proved that the In- 
dependent retarded integral in (2.8) is actually finite when B — > (the finite part 
is not followed by any pole). So is simply given by the value at B = of the 
retarded integral (see [8] and Section 4). But this is due to the special structure of 
the quadratic source Ntf" , and does not remain true at cubic and higher non-linear 
approximations (see for instance paper II). 

The first contribution solves (2.4a), but not the harmonic gauge condition 
(2.4b). The divergence of , say = d v u^ , is a priori different from zero. Using 
(2.6) we find 

w% = FP B=0 D^ 1 

The explicit factor B comes from the differentiation of r B in (2.8) (we denote 
rii = dir = x l /r). Owing to this factor the finite part in (2.10) is in fact a residue 
at B = 0, or coefficient of B~ x in the Laurent expansion. [The source term in (2.10) 
has a structure which is different from ', and unlike in (2.8) the integral admits 
in general a (simple) pole at B = 0.] The second contribution v% v in (2.7) is then 
defined in such a way as to compensate exactly the (a priori) non-zero divergence 
of , while being a homogeneous solution of (2.4a). This is possible because is 
a particular retarded solution of DwQ = (in the exterior region) . As such, it admits 
a unique multipolar decomposition in terms of four sets of STF tensors Al, Bl, Cl, 
Dl, namely 



B 



ro 



B 



rii 



(2.10) 



wl = Y,dL[r- 1 A L {t-r)] , (2.11a) 
' l =X>«- V~ X B L {t-r)] 

£>0 

+ {9l-i [r- 1 C lL - 1 {t - r)] + e iah d aL - 1 [r^D^^t - r)] } . (2.116) 



£>0 



These tensors can be computed straigthforwardly from the known expression (2.10), 
and the contribution v% u is defined in terms of these tensors. A particular definition 
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was proposed in the equations (4.13) of [20], where it was denoted q 2 v . Here we shall 
define this second contribution slightly differently, and accordingly we use the different 
notation v% v . The various components of v% v are given by 



v oo = _ r -ij A + da [ r -i (_ j Aa + jj Ca - SB a )] , (2.12a) 
vf = r- 1 (-JG + 3B 4 (1) ) - e iab d a [r^fD,] 

-J2 d L-i [r^AiL-t] , (2.126) 



£>2 



V 2 



+ d L -2 [r-\A^ L _ 2 + 3B^_ 2 - C ijL _ 2 ) 

- 2d aL - 2 [r- l e ah{t D j)hL _ 2 ] | . (2.12c) 

Like in (2.11) all the tensors are evaluated at the retarded time t — r. We note that 
the formulas (2.12) are non-instantaneous, as they depend on the moments Ml and 
Sl at any time less than t — r through the first and second time anti-derivatives 
of A, A a , C a , denoted e.g. by / A = f_Z A{t')dt' and ffC a = /*^ r / C a {t')dt' 
(see [20] for discussions). [To quadratic order the tensors Al, are given by 

some instantaneous functionals of the moments Ml and Sl-] The main property of 
is d v v 2 v = —W2, which is easily checked on the expressions (2.12). Furthermore 
Uv^ = 0, so the quadratic metric (2.7) is, indeed, a solution of both the wave equation 
(2.4a) and the gauge condition (2.4b). Note that the spatial trace v™ = 5^v l 2 J is 
especially simple, 



u« = Sr^B . (2.12d) 

The choice of definition (2.12) adopted here, which differs from the choice adopted in 
[20], is for convenience in future work. Of course we are free to adopt one definition 
or another because such a choice is equivalent to a choice of gauge. However, the 
definition (2.12) is slightly preferable to the definition proposed in [20] when we want 
to express the multipole moments Ml and Sl as integrals over the source. Thus we 
take the present opportunity to redefine the construction of the exterior metric using 
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this new definition. [Actually it can be checked that all intermediate and final results 
of this paper and paper II are independent of the choice of definition.] 

To the third (cubic) and higher non-linear iterations the construction of the 
external metric proceeds exactly along the same line, namely 

K v =<+< , (2-13) 

where the first term u£ u is the finite part of the retarded integral of the source to 
the nth post-Minkowskian order, ufa" = FY B=o^~fi[{r /tq) b N£ u {h\, h n -{j\, and 
where the second term vf^ v is defined from the divergence w% = d v u^ by the same 
formulas (2.11)-(2.12). (See [20] for the proof that the construction of the metric can 
be implemented to any post-Minkowskian order.) 



3. Structure of the quadratically non-linear field 

In this section we investigate the structure of the quadratic metric defined by 
(2.7)-(2.12). For simplicity we omit most of the numerical coefficients and indices, so 
as to focus our attention on the basic structure of the metric. A precise computation 
of the numerical coefficients is dealt with in Section 4. 

The structure of the linearized metric (2.3) is that of a sum of retarded multipolar 
waves, consisting of p spatial derivatives (say) acting on monopolar waves r~ 1 X(t — r), 

hi^aprt-r)] . (3.1) 

Our notation « refers to the structure of the expression. By expanding the derivatives 
dp (which act both on the pre-factor r _1 and on the retardation t — r) we get 

h x «^n Q r-''Z(*-r) , (3.2) 

where the powers of \jr are j > 1, and where we have expressed the angular 
dependence of each term using STF products of unit vectors uq = n < i i ni 2 ...rii q> = 
STF part of uq (see [31]). In practice, one may compute (3.2) from (3.1) by 

11 



decomposing dp on the basis of STF spatial derivatives 8q and using (A. 15) in 
Appendix A. After insertion of (3.2) into the quadratic source term one finds 

N 2 « J2 h Lr~ k F(t-r) , (3.3) 

k>2 

where the powers of 1/r start with k = 2 (as is clear from the fact that N 2 is quadratic 
in hi which is of order 1/r). The functions F are composed of sums of quadratic 
products of derivatives of the functions Z in (3.2). 

The main problem is to compute u 2 defined by (2.8). In view of the structure (3.3), 
it is a priori required to compute the finite part FP b=o of the retarded integral of any 
term hi,r~ k F{t — r) with multipolarity £ and radial dependence with k > 2. All the 
required formulas are listed in Appendix A (which summarizes results obtained mostly 
in previous works [20,8,9]). Notably, we know from Appendix A that the (finite part of 
the) retarded integral in the case k = 2 is irreducibly non-local or non-instantaneous. 
See (A.3)-(A.7) in Appendix A. When the power k satisfies 3 < k < £ + 2, where 
£ is the multipolarity (this excludes the monopolar case £ = 0), we know that the 
corresponding retarded integral is instantaneous, and given by (A.11)-(A.12). Finally, 
when k > £ + 3, we have again a non-local expression, given by (A.13)-(A.14), except 
in special combinations like (A. 16) for which the non-local integrals cancel out. 

The computation of u 2 = FPD^ 1 N 2 can be implemented by an algebraic computer 
program, following the successive steps (3.1)-(3.3) and applying on each terms 
composing (3.3) the formulas (A. 5), (A. 11) and (A.13)-(A.14). This is probably the 
most efficient way to obtain u 2 . However, in doing so we would discover that the only 
non-local integrals left in u 2 come from the source terms having a radial dependence 
with k = 2, in other words all the non-local integrals coming from terms having 
k > £ + 3 actually cancel out. Practically speaking the source terms having k > £ + 3 
turn out to combine into combinations such as (A. 16) yielding purely instantaneous 
contributions. This fact (proved in [8]) is special to the quadratic non-linearities, and 
does not stay true at higher orders, e.g. at the cubic order as seen in paper II. 

The proof of the latter assertion uses specifically the quadratic structure of the 
source N 2 . Instead of inserting in (2.5) the linearized metric in expanded form (3.2) 
and then working out all derivatives to arrive at (3.3), we keep the structure of the 
source as it basically is, i.e. a sum of quadratic products of multipolar waves, 
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N 2 ^^ j d P \r- 1 X{t-r) d Q r~ l Y(t-r)^ , (3.4) 

involving spatial multi-derivatives with P = %\...% p and Q = ji---j q and some functions 
of time X and Y. Then we perform on each term of (3.4) a sequence of operations by 
parts [i.e. diAdjB = di(AdjB) — AdidjB], by which the spatial derivatives acting on 
the wave in the left (say) are shifted in front and to the right. This leads to 

N 2 « J2 d p{ r ~ lx ( t ~ r)0fl[r -1 y(* - r)]} . (3.5) 

Only at this stage does one expand the space derivatives dp (inside the curly brackets), 
while leaving the derivatives dp in front un-expanded. The result reads 

N 2 ^ d P {n L r- k H(t-r)} , (3.6) 

2<k<£+2 

where the functions H are sums of products of X and time-derivatives of Y, and 
where we have projected the angular dependence of the terms inside the brackets on 
STF tensors ul- The point is that the radial dependence of the terms inside the 
brackets is related to the multipolarity £hy2<k<£ + 2. This can easily be seen 
from (3.5), as the expansion of the multipolar wave d^r^Y] is composed of a sum 
of terms 0L[r _1 y'] which have 1 < k < £ + 1 [see (A.l)], yielding 2 < k < £ + 2 
after multiplication by the factor r~ l on the left. The next operation is to single out 
the terms with pure radial dependence k = 2. All these terms can be obtained by 
applying dp on the terms inside the brackets having k = 2. In this way one generates 
besides all the terms k = 2 many other terms with k > 3, but the latter terms can be 
recombined into terms of the same form as in (3.6) (and thus having 3 < k < £ + 2). 
See [8] for the proof. Thus (3.6) can be rewritten as 

N 2 ^r- 2 Q(n,t-r)+ ^ 8 P {h L r- k F(t - r)} , (3.7) 

3<k<£+2 

where Q(n, t—r) denotes the coefficient of r~ 2 in the (finite) expansion of the quadratic 
source when r — > oo with t — r = const. The definition of Q(n, t — r) is 

*r=^(M-r) + 0(l) • (3.8) 
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Next, in anticipation of applying the finite part of the retarded integral, we multiply 
(3.7) by a factor r B , and introduce r B inside the brackets using again a series of 
operations by parts. In this way we get many new terms, but which all involve at 
least one factor B coming from the differentiation of r B during the latter operations. 
Thus 

r B N 2 « r B - 2 Q(n,t-r) + ^ 8 P {fi L r B ~ k F(t - r)} + 0{B) . (3.9) 

2,<k<t+2 

Applying the retarded integral on both sides of (3.9), commuting D^ 1 with dp, and 
taking the finite part, we are left with (the finite part of) retarded integrals of three 
types of terms: (i) the first term in (3.9) which has radial dependence r -2 , (ii) the 
terms in the brackets of (3.9) having radial dependence such that 3 < k < I + 2, 
and (iii) the terms O(B) which have the structure (3.3) with any radial dependence 
k but carry at least one factor B. In Case (ii), the retarded integrals are given by 
the instantaneous expressions (A. 11). In Case (iii), the retarded integrals are given 
by (A. 18) when the power of B is one, and are zero for higher powers. So in Case (iii) 
the retarded integrals are also instantaneous. Therefore we can state the following 
result [8]: the only non-local integrals in u 2 = FPD^A^ come from Case (i), namely 
from the source terms whose radial dependence is r -2 , and which are denoted by 
r~ 2 <5(n, t — r) in (3.8). These non-local integrals are given by (A.3)-(A.7). This 
result holds true only in the case of the quadratic non-linearity. In cubic and higher 
non-linearities, some hereditary integrals are generated by source terms with radial 
dependence such that k > £ + 3. 

Incidentally, note that the decomposition (3.9) of the source shows that the in- 
dependent retarded integral D B \ 1 [r B N 2 ] is finite at B = 0, i.e. does not involve any 
pole when B — > 0. Thus the finite part at B = is simply equal to the value of 
D~ 1 [r B N 2 ] at B = 0. This can be checked from the formulas (A. 3), (A.ll) and (A.18) 
in Appendix A, which are all finite at B = 0. Here again this is a peculiarity of the 
quadratic approximation. 

We are now in the position to write down the structure of the first contribution 
u 2 . From (3.9) and (A.3), (A.ll) and (A.18) we have 
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u 2 ^t 2 + ^2h L r k G{t - r) , 
fc>i 



(3.10) 



where the functions G(t — r) depend instantaneously on the multipole moments Ml 
and Sl (i-e. at the retarded time t — r only), and where the first term is non-local 
and given by 



On the other hand the second contribution v 2 , defined by (2.11)-(2.12), involves some 
time anti-derivatives of quadratic products of moments. We denote these time anti- 
derivatives by s 2 . Thus the quadratic-order metric can be written as 



where t 2 is the non-local integral (3.11), where s 2 are some anti-derivatives [given by 
(4.12) below in the case of the quadrupole-quadrupole interaction], and where we have 
many instantaneous terms. See (4.13) and Table 2 below for the complete expression 
of h 2 in the case of the quadrupole-quadrupole interaction. 



4. The quadrupole-quadrupole metric 

We specialize the previous investigation in the case of the interaction between two 
quadrupole moments M a b and M c d. Thus we keep in the linearized metric (2.3) only 
the terms corresponding to M ah , 



J1U _ 

2 _ 



a^[r- 2 Q^(n,t-r)] . 



(3.11) 




(3.12) 



hf = -2d ab (r^Mat) , 



(4.1a) 




(4.16) 



(4.1c) 



(Henceforth we use the same notation for the metric constructed out of the quadrupole 
moment as for the complete metric involving all multipolar contributions.) 
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Inserting (4.1) into the quadratic source (2.5), we can work out explicitly all the 
terms composing the source either in the all-expanded form (3.3) or in the more 
elaborate form (3.7). Notably, we find that the terms with radial dependence r~ 2 take 
the classic form of the stress-energy tensor of a massless field, 

Q^(n,t-r) = k»k v n(n t t-r) , (4.2) 
where /c M denotes the Minkowskian null vector /c M = (1, n l ), and where IT is given by 

n = n abcd M%> M% - An ah M^M^ + 2M%> M%> . (4.3) 

The quantity IT is proportional to the power (per unit of steradian) carried by the 
linearized waves. In the general case, Q^ v involves besides the quadrupole-quadrupole 
terms all the interacting terms between multipole moments with I > 2 and the mass 
monopole M. We introduce the STF multipole decomposition of IT, 

n(n,u) = ^n L U L {u) . (4.4) 

£>0 

From (4.3), the only non-zero multipolar coefficients are 



n = ^M^M^ , (4.5a) 

% = -f Mf>jt (4.56) 
U ijkl = . (4.5c) 

(n denotes the coefficient with multipolarity £ = or spherical average.) With this 
definition, and with the help of (A. 7) in Appendix A, we can write the non-local 
integral ti as 





P + OO 1 











-n(n',*-s). 4.6 

An s - rn.n' 

Some equivalent expressions follow from (A.3)-(A.5). For instance, we can write t% v 
in the explicit form 



r 



+ OO 

'M^Mifd^e} - AM^M^d%{A} + 2M ( $M ( $d^{2} 



(4.7a) 
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where the moments in the integrand are evaluated at the time t — s, where the multi- 
derivative operators mean for instance = d^d u didj with = (—d/ds,di), and 
where we use the special notation 

= ( 8 -r)Pln( 8 -r)-( 8 + r)Pln( 8 + r) 
plr 

[see (A. 4) in Appendix A]. 

Having the term , we undertake the computation of all the instantaneous terms 
x M c <i in [second terms in (3.10)]. The computation is straigthforward but 
tedious. As said above, when doing practical computations (notably by computer), 
the best method is the somewhat brute force method consisting of obtaining the 
source in the all-expanded form (3.3), and applying the (finite part of the) retarded 
integral on each term of (3.3) with k > 3 using the formulas (A. 11) and (A. 13). This 
is simpler than working out the source in the more elaborate form (3.9), and using 
the manifestly instantaneous formulas (A. 11) and (A. 18). The brute force method 
has also the advantage that one can check that all the non-local integrals but those 
coming from the r~ 2 term cancel out (as well as the associated logarithms of r). The 
term ui takes the form 



uf 



u 2 J 



& + E yu E \^n ahcd Mt k - m) M^ 

k=l m=0 ^ 

+ b^M^-^Mf + c k m Mf- k -^M^ } , (4.8a) 

fe=l m=0 ^ 

+ eth iab M£- k -^Mt } + ftnMt^M^ 

+ 9ina b M!- k - m) Mt } + htnMt k ~ m) M^ , (4-86) 

6 6—k s 

f ij , J_ J „k a . M (6-fe-m) M (m) 

k=l m=0 ^ 

+ ^n^Mf-^M^ + rt5 l3 h ahcd Mt k - m) 
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+ y k m n ab M£- k -^M$ + z k m M% k ~^M^ } , (4.8c) 

where all moments are evaluated at time t — r and where o^, b^, are purely 

numerical coefficients. Using the algebraic computer program Mathematica [32] we 
have obtained the numerical coefficients listed in Table 1. 

Next we follow the second part of the construction of the metric, and compute 
the divergence = d v u^ [see (2.10)-(2.12)]. To this end we need the divergence 
of the integral t^" , which is easily evaluated by noticing that dyt^ can be written, 
analogously to (2.10), as some residue at B = of a retarded integral (because of 
the explicit factor B), and furthermore that the radial dependence of the integrand is 
merely r~ s (because it comes from the differentiation of the source term r -2 ). From 
(A. 18), we know that when k = 3 the residue is non-zero only when the multipolarity 
is £ = 0. This yields immediately 



dj^ = FP B=0 n- 1 



D ( 



= -;/S fc " n(n -*- r) ' (4 - 9) 



The \x = component of (4.9) is proportional to the angular average of IT, already 
computed in (4.5a). One can check that the fx = i component is zero. Thus, 



d v t\ V = —r-'M^M^ , (4.10a) 
dvti? = . (4.106) 

Knowing (4.10) we can obtain by direct differentiation of the expressions (4.8), 
using the coefficients in Table 1. Again the computation is quite lengthy, but 
it provides us with an important check. Indeed, from (2.11) one must find that 
the divergence w% is a solution of the source-free d'Alembertian equation, namely 
□1^2 = 0. This test is very stringent, as a single erroneous coefficient in Table 1 would 
signify almost certainly its failure. Thus we determine all the tensors Al, Dl in 
(2.11). The non-zero ones are given by 



A = 



2 

75 



one ab ab Ar ab ab 
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An = -| {M { a %M J>a + M^Mfl) , (4.116) 



B = ~^M^M ab - -M%M% + , (4.11c) 

^ = M 3>a - \M^ a + S -Mi%M% a , ( 4.11 rf ) 

Di = e iab (-^Mg/Mtc - 2M£>m£> - ^M^M^ . (4.11e) 

By inserting the latter values into (2.12), we obtain the components of the second part 

v% v . In particular, we find some time anti-derivatives which define [see (3.12)] as 

A rt-r 

s oo 



Ir' 1 J duM$ Mf b \u) , (4.12a) 
J^duM^M^iu)^ , (4.126) 



82 - 5 

4 J = . (4.12c) 

The physical interpretation of these anti-derivatives is clear. Indeed the linearized 
metric (2.3) depends in particular on the mass monopole M and current dipole Si 
of the source (both are constant). Physically M and Si represent the total mass- 
energy and angular momentum of the system before the emission of gravitational 
radiation. Now the integral s® given by (4.12a) represents a small modification, due 
to the emission of radiation, of the initial mass M. This is clear from the comparison 
of (4.12a) and (2.3a), showing that there is exact agreement with the energy loss 
by radiation as given by the standard Einstein quadrupole formula. This result is 
originally due to Bonnor [10], and Bonnor and Rotenberg [4]. Similarly the integral 
s°* given by (4.12b) represents a modification of the total angular momentum in 
agreement with the quadrupole formula for the angular momentum loss. [There is no 
loss of total linear momentum at the level of the quadrupole-quadrupole interaction 
(one needs to consider also the mass octupole M a b c and/or the current quadrupole 

Sab)-] 

The quadratic metric can now be completed. We add up the two contributions 
u^ 1 [given by (4.8) and Table 1] and v^" [given by (4.11) and (2.12)]. The local terms 
are written in the same form as in (4.8). Thus, 

6 ^ 6— k s 

h% u = + + ^ — ^ < same expressions as in (4.8abc) but 



A— 1 ' m (I ^ 
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with coefficients a'£, , (4-13) 



where the non-local integrals and are given by (4.7) and (4.12), and where all 



the coefficients a'J^, ...,z'J^ are listed in Table 2. 



5. The quadrupole-quadrupole metric in the far zone 

We investigate the behaviour of the quadrupole-quadrupole field in the far 
zone, near future null infinity (i.e. at large distances when we recede from the source 
at the speed of light). The degrees of freedom of the radiation field, at leading order 
in the inverse of the distance, are contained in the transverse and tracefree (TT) 
projection of the spatial components ij of the metric (say <7^ T )- These are the so-called 
observable (or radiative) multipole moments, which are "measured" in an experiment 
located far away from the system. The TT projection gJ T to first order in 1/r reads 



tt T") J tt t/ I n 

9ij = - V Hab 2^ ^yiL-2U ijL -2 - j^n aL - 2 e ah{l V j)bL _ 2 ^ +V\^2 

(5.1) 

(with G = c = 1), where Ul and Vl denote the mass- type and current-type observable 
moments (both are functions of t — r) , and where the TT projection operator is 

Vijab(n) = (S ia - riin a )(5 jb - rijn b ) - ^(<% - ninj){5 ah - n a n b ) . (5.2) 

The ^-dependent coefficients in (5.1) are chosen so that Ul and Vl agree at the 
linearized order with the £th time derivatives of the moments Ml and Sl [compare 
with (2.3c)]. 

Let us consider first the non-local integral tQ" [see (4.6)-(4.7)]. As we know from 
(A. 8), the asymptotic expansion when r — > oo, t — r = const of the retarded integral 
of a source with radial dependence r~ 2 is composed of terms l/r n and lnr/r n . As 
such, behaves like In r/r when r — > oo, t — r = const. The logarithm is due 
to the deviation of the flat cones t — r = const in harmonic coordinates from the 
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true space-time null cones. The metric in harmonic coordinates is not of the normal 
Bondi-type at future null infinity [33,34]. Removal of the logarithm is done using 
radiative coordinates, so defined that the associated flat cones agree, asymptotically 
when r — > +00, with the true null cones (see e.g. [21]). This method, adopted in [9], 
permits to compute the observable moments in the non-local term ■ Here we follow 
another method, found by Thorne [15] and Wiseman and Will [16], which consists of 
applying first the TT projection operator (5.2) on . Because the TT projection 
kills any (linear) gauge term in the 1/r part of the metric, this method shortcuts the 
need of a transformation to radiative coordinates. However, one must be cautious in 
taking the limit r — > 00, t — r = const using (4.6). It is not allowed, for instance, to 
work out a leading 1/r term from the second expression in (4.6) because this term 
would involve a divergent integral (in accordance with the fact that the leading term 
is actually lnr/r). But, as pointed out in [15,16], the divergent parts of the integral 
cancel out after application of the TT projection, and at the end one recovers the 
correct result. Thus, we compute 

VgF 1 = -P«*V> f <h / f^n^f-,). (5.3) 

[Note that the TT projection as defined in (5.2) is purely algebraic. Strictly speaking 
it agrees with the true TT projection only when acting on the leading 1/r term.] With 
the multipole decomposition (4.4) we have 

«?r = -*W») g £' *. n,M / f • <s.4) 

We decompose the product of unit vectors n' a n' b n' L on the basis of STF tensors [31], 
we drop the terms having zero TT-projection, and we express the remaining terms 
using the Legendre function of the second kind Qi [see (A. 6c)]. Restoring the traces 
on the STF tensors, and droping further terms having zero TT-projection, we obtain 



(4 j ) TT = -^E c^ + ax^ + lU-i) ^- 2 / duU ^-^ 



x i(2e-l)Q e+2 [—) -2(2£ + 1)qJ — ) +(2£ + 3)Qi- 2 " "'" 





pt—r 


-2 






J — 00 


ft 


— u\ 


\ r J 



r 

(5.5) 
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The limit at future null infinity can now be applied. Indeed it suffices to insert the 
expansion of Qe(x) when x — > 1 + as given by (A. 8). As expected the limit is finite, 
because the terms \n(x — l)/2 in the expansions of the Qis cancel out. This yields 
[13,15,16,9] 



<*?) TT = -h^E ^/>"-w +° (£) ■ m> 

In the case of the quadrupole-quadrupole interaction we find [using (4.5)] 



(tl J ) TT = -V ijab [ du 

(where the brackets <> denote the STF projection). The non-local integral (5.7) 
represents the quadrupole-quadrupole contribution to the non-linear memory [14-16,9]. 

We now turn our attention to the instantaneous part of the metric. All the terms 
have been obtained in Section 4. We need only to apply the TT projection on the 
\jr part of as given by (4.13) and the coefficients listed in Table 2. After several 
transformations using STF techniques we obtain 



^ 1V1 c<a m b>c y§ Cd < ah cd > 



r(3) B/f (3) 



+ o 



lnr 



(5.7) 



(*? " ^) TT = \v ijab [- 2 -Mi%M b>c + + 4 -M^ a Mg\ 



+ n cd 



+ o 



10 

"I" Tldg^acg [ &ef<b 
1 



10 

'- \, 

10 ' 



-M^M d>f - \M^M^ f 



(5. 



(where the overbar on e means that the index e is to be excluded from the STF 
projection). 

With (5.7) and (5.8) we deduce the observable moments by comparison with 
(5.1). The quadrupole-quadrupole interaction contributes only to the observable mass 
quadrupole moment Uij, mass 2 4 -pole moment Uijki, and current octupole moment 
Vijk- We find 
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6U, 



7J-OC 



a<i J>a 



L a<i Ivl j> a 



a<i ]>a 



a<t J>a 



(5.9a) 



5Uijki — - 



t-r 



— M (3) M (2) 
5 IVI <ij IVI ki> > 



SVijk — £ a b<i 



— M^M k>b 



2 ja lv± k>b 



(5.96) 
(5.9c) 



Note that the non-local integrals are present in Uij and Uijki, but not in the current 
moment Vijk- 

Finally we add back in (5.9) the factor G and the powers of 1/c which are required 
in order to have the correct dimensionality. When this is done we find that SUij is 
of order 1/c 5 in the post-Newtonian expansion (2.5PN order), while both 5Uijki and 
SVijk are of order 1/c 3 or 1.5PN. This permits writing the complete expressions of 
Uij to 2.5PN order, and U^m, V%jk to 1.5PN order. In the case of U^ the reasoning 
has been done in [25], which shows that besides the quadrupole-quadrupole terms 
M ab x M c d there is an interaction between M ab and the (static) current dipole S c also 
at 2.5PN order, and there is the standard contribution of tails (computed in [9]) at 
1.5PN order. [In principle there is also an interaction between the mass octupole M abc 
and the mass dipole Ma at 2.5PN order, but we have chosen Ma = 0.] Thus, from 
(5.9a) and the equation (5.7) in [25], 



2GM /-*- r / c 



Uij(t-r/c) = M^f + - — - / a a 

G { 2 [ l - rlc , , r(3) (3 ) , , I 



/ 

J —c 



, . t — r/c — u\ 11 
M T + 12 



duM^M^u) + -Mi%M J>a 



(5.10) 

The tail integral involves a constant 11/12 computed in Appendix B of [9]. (See also 
[9] for the definition of the constant b.) The coefficient in front of the term M ab x S c is 
computed in Appendix B below. In a similar way, we find that U^ki and V^k involve 
the same types of multipole interactions, but that the terms M ab x M c d are of the 
same 1.5PN order as the tail terms. Thus, 
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U ijk i(t - r/c) = M\f kl + -S { 2M ! ' <lu 



, , t — r/c — u\ 59 

In '-, H 

2b J 30 



+ 



2 r 
5y_ 



oo 

* duM%M^{u)- 2 ^M%M kl> 



M ijli( u ) 



63 



102 



y\f ( 3 ) y\# ( 2 ) 
m <ij m kl> 



+ o 



1 



(5.11) 



, . t — r/c — u\ 5 

1,1 '—s— + 3 



+ £ ab<j 



1 



1 



rW»/(l) 



10 



<i m jk> } 



(5.12) 



The constants 59/30 and 4/3 are obtained from Appendix C in [24]. The coefficient 
of the term S < jM-^ > is computed in Appendix B below. 

The obtention in (5.10) of the observable quadrupole moment Uij up to the 2.5PN 
order yields the total power contained in the radiation field (or energy flux) complete 
up to the same 2.5PN order. Indeed it suffices to insert into (5.10) the intermediate 
quadrupole moment Mij determined in [25] as an explicit integral extending over the 
source at the 2.5PN order (the other needed moments being the octupole Mj jfc and 
current quadrupole Sij at 1PN order, and M^i and at Newtonian order, which 
are all known). Note that in order to obtain the wave form itself (and not only the 
energy flux it contains) one needs and Sij at the higher 2PN order. 



Appendix 1. Formulas to compute the quadratic non-linearities 

This appendix presents a unified compendium of formulas, many of them issued 
from previous works [20,8,9], which permit the computation of the quadratic non- 
linearities (involving any interaction between two multipole moments). The source of 
the quadratic non-linearities takes the structure (3.3), so we present the formulas for 
the finite part [as defined in (2.8)] of the retarded integral of any term fiLr~ k F(t — r) 
with multipolarity i and radial dependence r~ k where k is an integer > 2. [For 
practical computations it is convenient to use the source in expanded form (3.3) rather 
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than in the more precise form (3.9).] The problem was solved in [20] which obtained 
a basic formula for the B-dependent retarded integral 



□ 



R 



ro / r K 



r) 



>fe-3 



X 



r+oo 

I 



dsF(t- s)d L 



(2r ) B (B - k + 2)(B - k + 1) • • • (B - k - I + 2) 

( s _ r ^B-k+£+2 _ ^ _|_ r ^B-k+£+2 



(A.l) 



This formula is valid (by analytic continuation) for all values of B in the complex 
plane, except possibly at integer values of B where there is a simple pole. Note that 
to the STF product of unit vectors til in the left side corresponds a STF product of 
spatial derivatives 8l in the right side [31]. 

Our first case of interest is that of a source having k = 2. In this case (A.l) 
becomes 



□ 



R 



B 



x 



F(t-r) 



dsF(t-s)d L 



2(2r ) B B(B-l)---(B-e) 

s -r)B+£_ {s + r) B+£ 



(A.2) 



of which we compute the finite part in the Laurent expansion when B — > 0. We repeat 
briefly the reasoning of [8] : the coefficient in front of the integral admits a simple pole 
at B — 0, but at the same time the integral vanishes at B = thanks to the identity 
(A36) in [20] (see also (4.20a) in [8]). As a result the right side of (A.2) is finite at 
B = 0, in agreement with the fact that the retarded integral in its usual form (2.9) is 
convergent, with value 



□ 



R 



x di 



+oo 



dsF(t - s) 



(s — rY ln(s 



(s + rY hi(s + r) 



i\r 



(A.3) 



where we have removed the reference to taking the finite part at B = 0. Note that the 
length scale ro drops out in the result (this is thanks to (4.20a) in [8]). The formula 
(A.3) can be generalized to the case where the angular dependence is contained in any 
(non-tracefree) product k ai ...k ae of t Minkowskian null vectors k a = ( — l,n), 
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□ 



R 



kai ~: kae F(t-r) 



dsF(t - s) 



x d ai ... d, 



(s — r) £ ln(s — r) — (s + r) e ln(s + r) 



(A.4) 



where the space-time derivatives in the right side are defined by d a = (<9/<9s, <9;). A 
useful alternative form of (A. 3), proved e.g. in Appendix A of [9], reads 



□ 



R 



F(t-r) 



^ + °° dsF(t-s)Q e (^j , (A.5) 



where Qg(x) denotes the £th-order Legendre function of the second kind (with branch 
cut from — oo to 1, so x > 1). The Legendre function is given by 



Qt(x) 



' P,(y)-*- 
-i x-y 



\ / j=l J 



(A.6o) 
(A.66) 



where Pg is the Legendre polynomial (see e.g. [35]). Note that by combining (A. 6a) 
and the expansion of l/(x — n.n') in terms of Legendre polynomials (see e.g. (A26) 
in [20]), one has 



n L Qt(x) = f 



(A.6c) 



Atv x — n.n' 

where the angular integration dVt' is associated with the unit direction n n . Combining 
(A.5) and (A. 6c) we obtain another alternative formula, 



□ 



R 



1 



F(n,t-r) 



P + OO /• 



dfl' F(n',t- s) 



(A.7) 



An s — rn.n' 

This formula is valid for any function F(n, u) [not only for a function having a definite 
multipolarity I like in (A.5)]. It can also be recovered from the retarded integral in 
its usual form (2.9). 

The leading terms in the expansion when r — > oo (with t — r = const) follow from 
the expansion of the Legendre function when x — > 1 + . The expression (A. 6b) yields 
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Qt(x) = - 




(A.8a) 



and with (A. 5) this implies [9] 

>H = tf~ dAF(t - r - A) K^) + gj 




(A.86) 



This formula gives the leading term In r/r and the sub-dominant term 1/r. If we try 
to find the leading term using (A. 7) instead of (A. 5) [i.e. by changing the variable 
s = r + A in (A. 7) and expanding the integrand when r — > oo, t — r = const], we get 
formally a 1/r term but in factor of a divergent integral, as expected since the leading 
term is actually In r/r. 

In the case of a source term corresponding to k > 3, the relevant formula is 
obtained by performing k — 2 integrations by parts of (A.l). We get in this way 




□ 



R 



B 

r \ til 



ro 



B k-3 



= Or) 



i=0 



n L - 



B 

r \ n L 
ro 



F®(t-r) 



F( k -Q(t-r) 



(A.9) 



where the second term is a retarded integral of the type studied before, and where the 
coefficients are 



2 l (B- k + 2 + i)..(B- k + 3) 



(B-k + 2-£ + i)..(B -k + 2-£)(B-k + 3 + £ + i)..(B - k + 3 + i) ' 

(A.lOa) 

B(m = 2 k -iB(B-l)..(B-£) 

f { ] (B- k + 2)..{B- k-£ + 2)(B + £)..(B- k + £ + 3) ' 1 J 

The factors symbolized by dots decrease by steps of one unit from left to right. Before 
taking the finite part one must study the occurence of poles at B = in the coefficients 
(A. 10). Two cases must be distinguished. In the case 3 < k < £ + 2, none of the 
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denominators in (A. 10) vanish at B = 0. This implies that the second term in (A. 9) 
is zero at B = 0, owing to the explicit factor B in the numerator of /3(B). Thus we 
obtain (in the case 3 < k < I + 2) a finite result when B — > 0, which is local and 
independent of ro, 



□ 



Ft 



B 



2 k - 3 (k-3)l(e + 2-k)l 



(£ + k-2)\ 
>p (l + j) ! F( fc " 3 -J) (t - r) 



r j+i 



(A.ll) 



This formula can be checked by applying the d'Alembertian operator on both sides. 
When r — > oo, t — r = const, it gives 



□ 



F(t-r) 



2 k ~ 3 (k-3)l(£ + 2-k)K F( k ~ 3 \t-r) 

— TIL 



(£ + k-2)l 



(A.12) 



Next, in the case k > £ + 3, the denominators of both a(B) and /3(B) vanish at 
B = 0, so a(B) involves a (simple) pole, while (3(B) is finite (because the pole is 
compensated by the factor B in the numerator). In this case (k > £ + 3) the finite 
part reads 



FP B=0 



L ) B %m-r) 



^ F^(t-r) (-) k+e 2 k - 2 (k-3)l 
2_^o>in L t„ ~ n V 

i=0 



r k - 1 - 2 (k + £-2)\(k-£-3)\ 



x< (^ la ^)g t p (t - < ^(«-o 1+a -. 



The coefficients are given, when < i < k — £ — 4, by 



(A.13) 



(-2) l (A;-3)! (/c + £-3-z)!(/c-£-4-z)! 



(/e + £-2)!(/c-£-3)! 



(k-S-i)\ 



(A.Ua) 
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and, when k — £ — 3 < z < /c — 3, by 



2*(-)*+*(fc-3)! (k + £-S-i)\ 

r\ ■ == 

1 (k + £- 2)\{k - i - 3)! (A; - 3 - %)\{t - k + % + 3)! 

k-e-3 A £+3+i-k . A:-3 



X 



E 7- E 7- E 7+ E ]}■ (a. 146) 

j=l J 7=1 J j=k-2-i J j=k+£-2-i J } 



[One can express (A. 14) with the help of the Euler Y- function and its logarithmic 
derivative ip.] Note that (A. 13) has a genuine dependence on the length scale r 
through the logarithm of r/ro, which is in factor of a homogeneous solution, say 

A f CC-^ -i ^ V (£+j)! G«--»(t-r) 

The last term in (A. 13) involves a non-local integral known from (A.3)-(A.5) or (A. 7). 
There exists a special combination of retarded integrals (A. 13) which is purely local, 
finite at B = 0, logarithm-free and independent of ro- This combination, particularly 
useful in practical computations, reads 



□h 1 



B 

r 



) ^ ( 2 <* - 2) F< "'1 + K* - - 2) - V + l)]^P) 

=iI £M M ( ^-rn (A ;; 6o) 



where 



(-) fc 2 fc " 2 (A;-3)! 



7 



(fc + £ - l)!(fc 



2)! 



(fc + * - l)(fc - ^ - 2) - (2fc - 3)(Jfe - 2) 



(A.166) 



The formula is valid when k > £+2; when = £+2 we recover (A. 11). The second term 
in (A. 16a) is an homogeneous solution fixed by our particular way of integrating the 
wave equation. Thus the first term by itself is a solution of the required equation, but 
the homogeneous solution must absolutly be taken into account when doing practical 
computations with this method. 

To leading-order when r — > oo with t — r = const, (A. 13) reads 
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FP B= on^ 



— — F(t-r) 
r J r K 

• + CSO 



{k + t-2)\{k-i-S)\~ 

k-e-3 1 k+£-2 



X 



f + OO r / \ \ J 1 re-t-t — ^ -i 

/ dAF (-) (( _ r _ A) ln (M + E i + E i 



(A.17) 



As in the case k = 2 given by (A. 8b), the leading 1/r term is non-local. However, in 
contrast with (A. 8b), the expansion is free of logarithms of r (this is true to all orders 
in 1/r), and depends irreducibly on the length scale tq. 

Finally we give the expression of the pole part when B — > of the retarded integral 
in the case k > £ + 3. The result, easily obtained from (A.9)-(A.10), is 



□ 



-l 

R 



(-) fc 2 fc " 3 (A;-3)! 3 f F( k - £ ~ 3 \t - r) 



(k + £-2)\(k-£-3)\ 

(A.18) 

As there are only simple poles, the result is zero when the power of B is strictly larger 
than one. 



Appendix 2. The quadrupole- (current-) dipole interaction 

We start from the linearized metric (2.3) written for the mass quadrupole M ab and 
the (static) current dipole S a , 

/i?°=-2c> ab (r- 1 M a6 ), (B.la) 
hf = - 2e tab d a (r- 1 )S b + 2d a (r- 1 Mi 1) ), (B.16) 
h\ j =-2r- 1 i\4 2) , (B.lc) 

and replace it into (2.5), keeping only the terms M ab x S c . 

Once the source for this interaction is known in all-expanded form (3.3), we 
can apply on each of the terms the formulas of Appendix A. However, as we are only 
interested in the two coefficients of the terms M ab x S c in (5.10) and (5.12), it is 
sufficient to obtain the 1/r part of the metric when r — > oo, t — r = const. Thus 
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we proceed like in Appendix C of [24], which computed the 1/r part of h^" in the 
cases of quadrupole-monopole interactions. We recall the necessary formulas derived 
in [24], 



FP 



■ 1 d L (r- 1 F(t-r))]=( YflL 
t 1 



A 

2r 



2r 

+ o 



d\F^\t-r-\) 



lnr 



(B.2a) 



FP 



F^(t-r) 

= 2(£ + i) ( niflL ~ 6i < a * TlL - 1 >) ; + ° 



(B.26) 



FP B= on^ l^dijir-^dL (r^Fit-r)) 



--^SijSeo + (riij + Sij)h L 



2(£ + l)(£ + 2) 
- 2[5 i<ai n L -i > n j + dj <ai nL-i>ni] + 25 i<ae 5 jai _ 1 n L -2> 

+ o(X 



F(* +1 \t-r) 



(B.2c) 



We have added in (B.2c) the term — ^SijSio (where 5ij and 5eo are Kronecker symbols) 
with respect to the formula (C3) given in Appendix C of [24]. Indeed this term is 
mistakenly missing in the formula (C3) of [24] (but it does not change any of the 
results derived in [24]). 

The formulas (B.2) yield straightforwardly the 1/r term in the first part of 
the metric [see (2.8)]. Then we compute the (1/r term of the) divergence = d v u^ 
and deduce from (2.11)-(2.12) the second part v^" of the metric. We find that v^ v is 
non zero in the case M ab x S c , contrarily to the cases M x Ml and M x Sl studied 
in [24] for which the r^'s are zero. The result of the computation is 



h°2° = ~ 2 r ~ 1 e<*c n °dM$S c + 0(r" 2 ) , 



■^e abc n a M£'S c - -e iab n acd M^' S b 



r(4) 



+ 0(r" 2 ) , 



(B.3a) 
(B.36) 



h 



(4), 
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+ SijSabciiadM^ S c - 2e ah{i n ac M i ^ ) c S l 



(4) 

+ 2 £ ab(i n j)acd M c d bb 



+ 0(r" 2 ) . 



(B.3c) 



From (B.3c) we obtain the TT projection and then deduce the associated observable 
moments. Only the mass quadrupole and current octupole moments receive a 
contribution, 



SUij =-e ab<i M^ a S b , 



5Vi 



ijk 



- 2M%S k> 



(B.4o) 
(B.46) 



We thus obtain the coefficients quoted in (5.10) and (5.12) of the text. 
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Table 1: The numerical values of coefficients entering the expression of u^" as defined in (4.8). 
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